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We explore the microwave radiation emitted from a biased double quantum dot due to the inelastic 
tunneling of single charges. Radiation is detected over a broad range of detuning configurations 
between the dot energy levels with pronounced maxima occurring in resonance with a capacitively 
coupled transmission line resonator. The power emitted for forward and reverse resonant detuning 
is found to be in good agreement with a rate equation model, which considers the hybridization of 
the individual dot charge states. 


The electronic properties of semiconductor nanostruc¬ 
tures have been widely studied using transport measure¬ 
ments [1]. More recently, the use of radio and microwave 
frequency measurement techniques has enabled and stim¬ 
ulated a new generation of experiments [2-4] mostly with 
superconducting qubits. In similar device geometries the 
coupling of double quantum dots in carbon nanotubes 
[5], gate-defined GaAs heterostructures [1, 7], InAs nano¬ 
wires [8], and graphene [9] to GHz frequency coplanar 
waveguide resonators has been explored. Earlier exper¬ 
iments have investigated photon emission and lasing ef¬ 
fects in superconducting circuits [10, 11]. More recently 
experiments studied photon emission from biased dou¬ 
ble quantum dots. Gain and micromaser action have 
been predicted [12] and observed by pumping a single 
microwave resonator mode [13, 14]. 

Here, we experimentally explore inelastic tunneling in 
a semiconductor double quantum dot capacitively cou¬ 
pled to a transmission line resonator by detecting the 
microwave radiation emitted in the process. The de¬ 
tection of the weak microwave signals is facilitated by 
the use of near quantum limited parametric amplifiers 
[15, 16]. Previous experiments on quantum dots coupled 
to microwave resonators [1, 5, 7-9] detected mostly po¬ 
larizability allowing to extract charge stability diagrams. 
Here, we demonstrate that the level separation of hy¬ 
bridized states in a double quantum dot, if on resonance 
with the microwave resonator, can be investigated with 
high precision by directly detecting inelastic transitions. 
In particular, the tunneling rates, which are typically dif¬ 
ficult to measure in pure dc transport experiments, are 
directly reflected in the power of the respective emitted 
microwave radiation. Future options include the possibil¬ 
ity to characterize the classical and quantum properties 
of radiation emission from semiconductor nanostructures 
at microwave frequencies using correlation function mea¬ 
surements [17] or state tomography [18]. 

The hybrid device explored in our experiments is real¬ 
ized by capacitively coupling a gate-defined double quan¬ 
tum dot to an on-chip superconducting coplanar waveg¬ 
uide resonator [Fig. 1]. The gate structure depicted in 
Fig. 1(b) is patterned on top of a GaAs/AlGaAs het¬ 
erostructure with a two-dimensional electron gas (2DEG) 
90 nm below the surface. The quantum dots are formed 



FIG. 1. (a) Energy level structure of a double quantum dot 
with levels eu and e_B. The dots are tunnel coupled to source 
(S) and drain (D) with chemical potentials ps and po at rate 
Tl and Tr. The interdot coupling rate is given by t/h. We 
consider the creation of photons in a cavity of resonance fre¬ 
quency Vo in dependence on the interdot detuning S = er, — cr. 
In the depicted situation (<5 > 0) the detuning is in forward 
direction with respect to the source-drain bias, (b) Scanning 
electron micrograph of the quantum dot structure: center gate 
(CG), left and right plunger gate (LPG, RPG), quantum point 
contact (QPC). The LPG shaded in orange capacitively cou¬ 
ples the left dot to the resonator (Glpg). (c) Schematic of 
sample and measurement setup. A coherent microwave signal 
(vr) is applied to the resonator. The output signal passes a 
circulator and is amplified by a Josephson parametric ampli¬ 
fier (JPA) and a high electron mobility transistor (HEMT) 
amplifier. It is mixed with a local oscillator signal (r'Lo) in a 
heterodyne detection scheme and processed with a field pro¬ 
grammable gate array (FPGA). 


by negatively biasing the top gates. The left and right 
plunger gates (LPG, RPG) control the electrochemical 
potentials (eL,eR) of the electrons in the respective dots 
and thereby the detuning energy 5 = ep — cr between the 
dot levels. The center gate is used to adjust the interdot 
tunnel coupling energy t. The tunnel coupling leads to 
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FIG. 2. (a) Change in phase A(f> of a constant RF measure¬ 
ment tone applied to the resonator as a function of left and 
right plunger gate voltages in the vicinity of the triple points. 
The borders between stable charge configurations with n elec¬ 
trons in the left and m electrons in the right dot are indicated 
by dashed lines. The arrow indicates the axis along which the 
detuning S between the dots is varied. A fluctuating back¬ 
ground constant in LPG voltage is subtracted, (b) Associ¬ 
ated shift in the resonance frequency Ai/ = v — vq against 
detuning 5. The solid line is a Master equation simulation of 
a Jaynes-Cummings model from which the coupling strength 
g/2n « 13 MHz, decoherence rate = 250 MHz and tun¬ 
nel coupling 2t/h ~ 4.4 GHz are extracted [1, 7]. 


a hybridization of the left and right dot states. The re¬ 
sulting bonding and antibonding states form a two-level 
system with an energy separation = \J ( 2 t )2 -|- P-^ 
where h is Planck’s constant and denotes the transi¬ 
tion frequency of the charge qubit [19]. 

The 200 nm thick Al coplanar waveguide resonator has 
a bare resonance frequency of uq = 6.852 GHz (28 fieV / h) 
and a loaded quality factor of 2058. The center con¬ 
ductor of the resonator extends to the left dot form¬ 
ing the plunger gate shaded in orange, which mediates 
dipole coupling between double quantum dot and res¬ 
onator [Fig. 1(c)] [1, 7]. The sample is similar to the 
devices presented in Refs. 1 and 20. 

We first characterize the properties of the coupled 
system in microwave transmission measurements as de¬ 
scribed in Ref. 1. This yields all parameter values rel¬ 
evant for the analysis of the photon emission data pre¬ 
sented subsequently. A microwave tone of constant fre¬ 
quency Vr set to the bare cavity frequency i>o is trans¬ 
mitted through the resonator. The output field is first 
amplified by a Josephson parametric amplifier [21] pro¬ 
viding quantum limited amplification [15, 16]. A high 
electron mobility transistor amplifier provides a second 
amplification stage before the signal is mixed with a 
local oscillator microwave tone. The signal is subse¬ 
quently digitally down-converted and processed with a 
field-programmable gate array, giving access to its am¬ 
plitude A and phase (j) [Fig. 1(c)] [ 2 ]. 

Charge delocalization between the dots leads to a dis¬ 
persive frequency shift of the cavity seen as a change in 


phase A(j) of the transmitted tone as plotted in Fig. 2(a). 
The charge stability regions are indicated together with 
the i5-axis along which the detuning between the dots is 
varied [22]. In the presented experiments, the number of 
electrons in each dot is on the order of 10 as determined 
by quantum point contact charge detection [23, 24]. 

The dispersive frequency shift as a function of 6 is 
measured by recording full transmission spectra. At de¬ 
tuning energies in the vicinity of the qubit-cavity res¬ 
onance condition the resonator frequency is strongly 
shifted [Fig. 2(b)]. We calculate the frequency shift nu¬ 
merically using a Jaynes-Cummings model [25] to extract 
the resonator-dot coupling strength g, the interdot tun¬ 
nel coupling energy t and the decoherence rate 7 . To 
analyze the emission experiments performed at different 
tunnel coupling energies t discussed in the following, we 
have used a constant value of g/2TT « 11 MHz, see sup¬ 
plemental material [36]. 

With this device, decoherence rates as low as 250 ± 
50 MHz were observed, which is more than an order of 
magnitude lower than for previous samples investigated 
in our group [ 1 , 20 ] and comparable to the lowest rates 
observed in other experiments [7, 26, 27]. We attribute 
the improvement to a combination of optimized filter¬ 
ing and a very stable 2DEG. It also manifests itself in 
a reduced electron temperature of ^ 60 mK as extracted 
from conductance resonances. Note that the improved 
coherence facilitates the emission experiments presented 
below. 

To study photon emission from the double dot, we ap¬ 
ply a bias of Vsd = —200 fiV (48 GHz) between source 
and drain. The bias is chosen to be smaller than both 
the charging energies and the single-particle level spacing 
to ensure single-level transport. The charging energy of 
each dot is roughly 1 meV (« 200 GHz ■ h), as extracted 
from the stability diagram, and at Vsd = —200 /rV trans¬ 
port measurements show no excited states in the bias 
window. The tunnel coupling to the leads, adjusted to 
Fl ~ Fr ^ 1 GHz, controls the current through the 
device. 

We measure the power spectral density (PSD) of pho¬ 
tons emitted from the cavity [2, 36] with the double quan¬ 
tum dot being the only source of radiation. For detuning 
values \6/h\ larger than the resonance frequency, we find 
a low photon signal of less than 0.01 Hz“^s“^, see purple 
curve in Fig. 3(a). Glose to resonance, the photon signal 
is significantly increased to more than 0.06 Hz“^s“^ (red 
curve) corresponding to an average number of 0.015 pho¬ 
tons in the resonator extracted from the integrated PSD. 
Both data sets fit well to a Lorentzian lineshape (solid 
lines) taking into account the frequency dependence of 
the parametric amplifier gain. The linewidths extracted 
from all PSD measurements are identical to the resonator 
linewidth k/2tt = 3.3 MHz. 

We measure the PSD vs. detuning 6 around the lower 
triple point shown in Fig. 2(a). Integrating Lorentzian 
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FIG. 3. (a) Power spectral density (PSD) of microwave emission measured at the two detnning energies 5 indicated by dashed 
circles in (b). The extracted linewidth corresponds to the resonator linewidth k. (b) Plot of the source drain cnrrent (blue, 
scale on the right hand side) and the emission power in units of photons emitted into the cavity per microsecond (orange) 
against detuning 5. Each point in emission power is obtained by integrating the respective PSD. (c) Photon emission power 
measured for the indicated interdot tunnel rates 2t/h. The background proportional to the current is subtracted and the 
emission resonances are fitted using a sum of two Gaussian lineshapes to extract the values of resonant detuning and power 
(sum: solid orange line; individual lineshapes: dashed gray lines). 


fits to the data yields the number of photons emitted 
from the cavity per unit of time, P. The power P is 
plotted in Fig. 3(b) as a function of detuning 5 (orange) 
and compared to the simultaneously measured source- 
drain current (blue). We observe a pronounced resonance 
in the measured power P at positive detuning (-l-(5o) as 
well as a lower peak at negative detuning (—(5o). These 
maxima occur when the energy released in the interdot 
transition corresponds to the resonator frequency {vg = 
vq) yielding the resonant detuning energies 

±(5o = ±\/(hr'o)^ - (2t)^ . (1) 

The background emission power away from the reso¬ 
nances is proportional to the current. We associate it 
with combined photon/phonon processes and tunneling 
processes between the left or right dot and the continuum 
of states in the leads. The proportionality factor is the 
same for all measurements and yields a background rate 
of roughly 1.3 x 10“^ photons emitted from the cavity per 
transported electron. This shows that competing relax¬ 
ation channels for the qubit, such as phonon emission, 
are relevant in inelastic processes even in the presence 
of the resonator [13, 29-31]. Irrespective of the details 


of the non-radiative relaxation mechanism we expect the 
rate of photon emission into the resonator to be limited 
by the ratio of the qubit dephasing rate to the cavity 
linewidth ~ 0.01. 

We investigate emission for interdot tunnel coupling 
rates up to 2t/h = 7.4 GHz. Representative examples 
are shown in Fig. 3(c). We find that the emission power 
at —(5o increases with t, while the separation between the 
resonances, indicated by vertical dashed lines, decreases. 
To analyze the resonances in emission in detail we sub¬ 
tract the background signal proportional to the current. 
A sum of two Gaussian lineshapes is fitted to the res¬ 
onances to analyze the maximum power and resonant 
detuning [Fig. 3(c)]. We define the zero detuning to be 
centered between the resonances. For details of the PSD 
measurement and the extraction of the tunnel coupling, 
refer to the supplemental material [36]. 

The extracted resonant detuning l(5ol//i approaches 
zero as 2t approaches hi/Q [Fig. 4(a)]. For 2t > h^o, 
we find a single resonance [Fig. 3(c)], which decreases 
in power and eventually vanishes when t is further in¬ 
creased. For 2t/h < A GHz the emission power at —(5o is 
below the noise level of our detection system for the cho- 
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FIG. 4. (a) Resonant detuning |5o| and (b) power ratio Pr/Pf 
as a function of tunnel coupling 2t. The data are compared 
to expressions (1) and (5) (solid curves). Error bars in |5o| 
and P-cjPi are determined by fitting error of the Gaussian 
lineshapes shown in Fig. 3(c), error bars in 2t result from 
the Jaynes-Gummings simulations of transmission measure¬ 
ments. (c) Level diagrams depicting the wavefunctions of 
the hybridized dot states for positive and negative detuning. 
The dashed black lines represent the dot energy levels in the 
left-right {\L),\R)) basis, the solid red lines represent the hy¬ 
bridized states and the shaded gray areas the wavefunctions 
of the bonding and antibonding state. 

sen integration time, so that no points can be acquired 
in this region. The experimental data are well described 
by Eq. (1) and deviations lie within the error margins 
[Fig. 4(a)] where all parameters are directly determined 
from the experiment. 

For the negative resonance condition — (5o, the detun¬ 
ing of the electronic states is reversed with respect to the 
source-drain bias [Fig. 4(c)]. In this configuration, trans¬ 
port is expected to be blocked if the hybridization of the 
occupation states \L),\R) of the left and right dots is not 
considered. The hybridized states are given by 

Je) = cos(6»/2)lL) -b sin{9/2)\R), (2) 

\g) = sin(6»/2)lL) - cos(6»/2)li?), (3) 

where 0 is given by tan0 = 2t/6. Considering the anti¬ 


bonding wavefunction je), we find the electron in the left 
dot with probability a = cos^( 0 / 2 ), which is a measure 
of the hybridization of the qubit states [Fig. 4(c)]. Even 
for reverse detuning, electrons can therefore tunnel into 
the excited state leading to photon emission at —do- 

We find that emission power for forward detuning Pf 
is roughly independent of t while photon emission for re¬ 
verse detuning P,- increases with t [Fig. 3(c)]. The power 
ratio Pr/Pf approaches unity as 2t approaches and 
the peaks merge [Fig. 4(b)]. 

We model the behavior of the emitted power with in¬ 
creasing t using a master equation approach. We describe 
the system by three possible states: qubit in the excited 
state je) (excess electron in antibonding state), qubit in 
ground state \g) (excess electron in bonding state) and 
qubit in state jO), which we define as the state with no 
excess electrons in the dots. This describes the system 
around the lower triple point. For the upper triple point, 
the third qubit state corresponds to two extra electrons 
in the dots leading to the same result. We arrive at the 
following rate equation for the occupation probabilities 
of the three qubit states [32]: 

Po= -rnPo + Oi^RPg + (l-a)rRPe 

Pg= (l-a)rLPo - aT^pg (4) 

Pe= aVi^po - (l-Q!)rRPe- 

Here, we neglect the relaxation rate 71 from the excited 
to the ground state. This is justified because we estimate 
7 i « 100 MHz [33], which is significantly smaller than 
the tunneling rates to the leads Tr « Tr. The same 
holds for the resonator-dot coupling g « 11 MHz. We 
confirmed that including 71 in the rate equation does 
not significantly affect the power ratio even for 71 ^ 
100 MHz. 

Because the resonant emission signal is generated in 
transitions from the excited state we expect the emission 
power to be proportional to Pe{5), which is obtained using 
the steady-state condition p = 0. This yields the emission 
power ratio 

_ Pej-Sp) _ {hv^Y - 2t^ - 

Pf Pe{+5p) {hvpY - 2P -b huo5o ' 

The experimental results are well approximated by this 
model [black solid line in Fig. 4(b)]. There is, however, 
a systematic deviation, such that the power ratios in all 
experimental data are slightly higher than the model pre¬ 
dicts. We can account for this deviation by considering 
energy dependent tunnel rates Tl and Tr, see Ref. 36. 
Note that, as for the resonant detuning [Fig. 4(a)], there 
is no free fitting parameter. 

In conclusion, we investigated photon emission from 
a biased double quantum dot dipole coupled to a mi¬ 
crowave cavity. We clearly distinguish between back¬ 
ground emission and resonant photon emission of the 
qubit generated in inelastic interdot transitions. We 
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found photon emission for both forward and reverse in¬ 
terdot detuning with respect to source-drain bias reflect¬ 
ing the distribution of wavefunctions in the hybridized 
electronic states. We modeled the results using a mas¬ 
ter equation approach for the double quantum dot and 
found good agreement with the experiment. In future ex¬ 
periments, microwave correlation function measurements 
could be used to investigate photon statistics [17, 34]. For 
the low photon rates we detect from single-electron inter¬ 
dot transitions in our experiments photon antibunching 
is predicted [35]. 
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MEASUREMENT DETAILS 

The emission experiment is performed at finite source- 
drain bias where the triple points of a charge stability 
diagram expand into finite bias triangles [Fig. SI]. We 
measure the emitted power at the lower triangle and vary 
the interdot detuning by tuning the left and right plunger 
gate voltages along the J-axis indicated by the dashed 
purple line [Fig. 3(b)]. We additionally study emitted 
power vs. tunnel coupling [Fig. 3(c)] by tuning the center 
gate [Fig. 1(b)] voltage. 
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Extraction of tunnel coupling 

For each center gate voltage, the tunnel coupling en¬ 
ergy is extracted individually from microwave transmis¬ 
sion measurements [1]. Without source-drain bias, trans¬ 
mission measurements show a phase shift Acf of the out¬ 
put microwave signal in the region of the charge degen¬ 
eracy line connecting the triple points [Fig. S2(a)]. At 
finite bias, the phase shift extends further along the base 
lines of the finite bias triangles. We aim at extracting 
the interdot tunnel coupling t of the biased system at the 
lower finite bias triangle corresponding to the conditions 
of the emission measurements [purple line in Fig. SI and 
Fig. S2(b)]. However, we observe that the charge trans¬ 
port occuring within the finite bias triangles changes the 
phase and frequency response of the resonator compared 
to the unbiased system. We therefore extract the tunnel 
coupling using the Jaynes-Cummings description from 
measurements without source-drain bias where there is 
no transport through the system. 

The tunnel coupling between the dots is not constant 
along the charge degeneracy line but decreases for more 
negative plunger gate voltages which is observed in the 
phase shift [Fig. S2(a)]. This is due to the effect of the 
cross-coupling of the plunger gates on the central tunnel 
barrier and thereby on t. We can quantify the change of 
tunnel coupling by measuring full transmission spectra 
vs. detuning for multiple cuts through the charge degen¬ 
eracy line [Fig. S2(a) and (c)] and hnd that t decreases 
linearly with plunger gate voltage. 

In the unbiased case the cavity frequency shift disap¬ 
pears at the triple points, so that the Jaynes-Cummings 
Hamiltonian cannot model the frequency shift Aiz and 
t has to be extracted between the triple points [dashed 


FIG. SI. Finite bias triangles at Vsd = —200 /rV showing 
the source-drain current / vs. left and right plunger gate 
(LPG, RPG) voltage. The dashed purple line marks where 
the emission data is taken. 

green line. Fig. S2(a)]. For the emission measurements 
presented in Fig. 3(c), the bare interdot tunnel coupling is 
extracted from the frequency shifts Ah' shown in Fig. S3 
and the values are given in the dashed green box. We 
extrapolate the linear relation between t and the plunger 
gate voltages to obtain t at the triple point. 

Additionally, we need to take the effect of the bias 
on the tunnel coupling into account. Applying a bias of 
Vsd = —200 fiV shifts the charge degeneracy line slightly 
in the LPG/RPG plane leading to a change of « 50 MHz 
in 2t/h for points on the charge degeneracy line [cf. or¬ 
ange curves in Fig. S2(c) and (d)]. With this correction 
we find the shifted 2t/h for the parameter settings of 
the emission experiment, which are listed in the dashed 
purple box in Fig. S3. 

The simulation based on the Jaynes-Gummings model 
yields 2t/h up to a systematic uncertainty between 50 
and 200 MHz, depending on the tunnel coupling. An 
additional error is introduced by the linear extrapolation 
described above. Finally, the correction for the bias shift 
also entails a small additional uncertainty. The error 
bars plotted in Fig. 4(a) and (b) include these three error 
sources. 

Power spectral density 

The held quadratures I and Q of the emitted mi¬ 
crowave radiation are measured in a heterodyne detec¬ 
tion scheme. Fast signal processing with a held pro- 
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(a) bias off (b) bias on 



(c) (d) 




FIG. S2. Measured phase shift A</) at Vsd = 0 (a) and 
VsD = —200 pV (b). The dashed black line serves as a ref¬ 
erence to illustrate the shift of the charge degeneracy line in 
the LPG/RPG plane when the bias is turned on. The dashed 
purple line indicates the gate voltage configuration where the 
emission data is acquired, while the dashed green line indi¬ 
cates where the bare tunnel coupling is extracted, i.e. where 
the curves shown in Fig. S3 are taken, (c), (d) Shift of the 
resonant frequency Au (colored dots) measured along the cuts 
indicated in (a) and (b). Solid lines are results of numerical 
simulations. In the left panel we observe that the tunnel rate 
2t/h (insets) increases with plunger gate voltage. We com¬ 
pare the middle (orange) cut between bias off and on and see 
that 2f changes by ~ 60 MHz in this case. 


grammable gate array (FPGA) gives access to correla¬ 
tions of the emitted photons. We compute the power 
spectral density by taking the Fourier transform of the 
first-order auto-correlation function of the complex signal 
amplitude S = I + iQ = Ae'‘^ [2]. 

We need to resolve a small photon signal on the order of 
10 “^ generated by the biased double quantum 

dot. To distinguish this signal from the noise, we subtract 
zero-bias reference data. For this purpose, we alternate 
between measurements with the bias turned on and off 
in this case taking 40 repetitions. After averaging the 
measurement (bias on) and the reference (bias off) data, 
we obtain the PSDs S'on and 5'off- Subtracting the two 
yields the power spectral density of the signal modes in 
which we are interested: PSD = Sa = Son — Sqs- Details 
about the measurement technique and signal processing 



I 2t/h at ; 
I emission; 
I cond. = 1 

I 4.0 GHz I 
i 4.7 GHz i 
I 5.9 GHz I 

; 6.3 GHz ; 

; 7.0 GHz ; 


FIG. S3. Measured frequency shift Av (colored dots) and nu¬ 
merical simulation (solid lines) to extract t. The bare values 
of 2t/h extracted from the curves are indicated in the dashed 
green box on the left. The values given in the dashed purple 
box on the right take into account the shifts described in the 
text and are valid at the conditions of the emission experi¬ 
ments presented in the main text. The other parameters of 
the simulation are the coupling strength g/2'K = 11 MHz and 
decoherence 7/271 = [500, 600, 800,1000,1000] MHz (from top 
to bottom). 


can be found in Ref. 3. 

The PSD plotted in Fig. 3(a) is scaled to represent 
the number of photons emitted from the resonator per 
Hz and per second. To obtain this scaling , the number 
of noise photons at the input of the parametric ampli¬ 
fier was determined using the AC Stark shift observed 
in a superconducting qubit measurement with a similar 
detection chain in the same cryostat. 


RATE EQUATION CALCULATIONS 

In the model presented in the main text we neglect 
the decay rate 71 and assume energy independent tunnel 
barriers. This yields the rate equation p = Mp for the 
probabilities to find the system in one of its three states, 
where p = {po,Pg,PeV and 


/ -Fl aFu (l-a)rR \ 

M= (l-a)rL -oFr 0 

V aFu 0 -(l-a)rR/ 

We calculatepe(^) using the steady-state condition;? = 





























0 and po + Pg + Pe = 1 and find 


(a) 


Pe{S) 


Fl {2t^ + S {S + VWTP)) 

2t^TR + 2ri^)+2Ti^6^ 


Pe(±l5o) 


Tl {{hiyo)'^ - 2t^ ± hvoSo) 
2t^{T^-2Ti^)+2Ti,{hvo?' 


( 6 ) 


This expression describes the t-dependence of the power 
emitted at each resonance condition ±(5o and is indicated 
with dashed lines in Fig. S4(a). We assume Fl = Fr 
and scale the curves to agree with the measured emis¬ 
sion power. Experimentally, we observe that the emitted 
power at —Jq (blue) increases with t, while the emission 
power at -|-5o (red) hardly depends on t, which clearly de¬ 
viates from the behavior predicted by Eq. (6). A possible 
^-dependence of the decay rate from excited to ground 
state related to photon emission into the resonator is not 
taken into account. Note that the peaks are merged for 
the last two data points, which are therefore omitted in 
the figure. 

This deviation is also apparent when we consider the 
power ratio 

^ _ Pe(-<^o) _ - 2.t^ - hypSp 

Pi Pe(+i5o) - 2.P + hi/oSo ’ 

which is plotted in black in Fig. S4(b). No assumptions 
for the tunnel rates to the leads (Fl, Fr) and no scaling 
factor enter this expression. 



FIG. S4. (a) Emitted power at —do (blue) and -|-do (red). The 
dashed lines are the theoretical expressions given by Eq. (6), 
where we assumed Fl = Fr and a scaling factor. The solid 
lines are plots of pe(—do) (blue) and Pe(-l-do) (red) assuming 
F = r® = Fr = F^ and Fr = cF with c = 2.6. (b) Power ra¬ 
tio Pr / Pi- The solid black line shows the theoretical ratio ob¬ 
tained with energy independent tunnel barriers [cf. Fig 4(b)], 
while the solid orange line was obtained from the generalized 
model assuming energy dependent tunnel barriers. 


Energy dependent tunnel rates to the leads 

We can generalize the model by allowing an energy 
dependence of the tunnel rates Fr and Fr, i.e. the tunnel 
rates into and out of the excited state F£, Fr differ from 
those into and out of the ground state F®, Fr [4]. The 
energy dependence can arise from changes in the tunnel 
coupling or the density of states in the leads. We find 
the generalized transition matrix: 

/-(1 - a)F[ - aF£ oF® (1 - a)Ff, \ 

M = (1 - a)Ff -aF® 0 

V ar£ 0 -(l-a)F^y 

The steady state condition p = Mp = 0 now 

yields a modified excited state occupation probability 
Pe(< 5 , Fr, F£, Fr, Fr). A first assumption could be that 
the left and right barriers are both energy dependent in 
the same way, i.e. F£ = cF® and Fr = cFr with some 
factor c typically > 1. This, however, does not influence 
the theoretically predicted power ratio. 

We can correct for the deviation in Fig. S4(b) if we 
assume that only the tunnel rate from the left lead in¬ 
creases with energy: F£ > Fr. More specifically, we 


assume F = F® = Fr = Fr and Fr = cF. We fit the 
data to the modified power ratio yielding c = 2.6. The 
order of magnitude of this energy dependence seems re¬ 
alistic considering the findings of other transport exper¬ 
iments [5, 6]. We assume that the origin for the energy 
dependence lies in a modulated density of states due to 
gate geometry and disorder [6]. In the measured device 
the left and right side gate influence the double dot con¬ 
finement potential very differently, which is not expected 
from the device geometry. This could be the reason for 
the left/right asymmetry of the tunnel rates to the leads. 

The resulting modified expressions Pe(±<5o) for the 
emitted power and the power ratio Pr/Pf are shown in 
Fig. S4 and the agreement with the experimental data 
is considerably improved. Simply assuming Fr > Tr 
does not influence the power ratio, the asymmetric en¬ 
ergy dependence is therefore necessary to account for the 
deviation. 


Resonance linewidth 

Further information about the emission process can 
be gained from analyzing the linewidth of the emission 
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FIG. S5. Full width at half maximum (FWHM) of the emis¬ 
sion resonances at —5res (blue) and H-dres (red). 


resonances. The Gaussian peaks fitted to the data in 
Fig. 3.(c) have a full width at half maximum (FWHM) 
between 1 and 5 GHz increasing with t [Fig. S5]. 

This linewidth can be set in relation to the energy 
broadening of the electronic dot states. We assume that 
the resonances of the electronic states are Lorentzian 
with a FWHM r/2. The emission signal is described by 
a convolution of the two states involved in the transition, 
which yields a Lorentzian in Vg with linewidth F: 


pi^,) = m 


r /2 


The emission amplitude 7i(5) depends on the detuning 
between the dot states and is in our model proportional 
to the excited state occupation probability: H(i5) (x PeiS). 
The resonator linewidth k/2tt = 3.3 MHz is much smaller 
than F/2, so that the resonator response can be approxi¬ 
mated as a (5-function and neglected in this consideration. 

We are interested in the emission resonances at con¬ 
stant t as a function of detuning S: 


P{S/h) = H((5) 


r /2 

{^{2ty + 6yh-iyo)^ + {T/2)P 


(7) 


This is not Lorentzian in S but yields two peaks sym¬ 
metric in 6, explaining the two resonances we observe in 
emission. 


By fitting Eq. (7) to the measured emission resonances, 
we hnd F = 1.5 ± 0.4 GHz. This corresponds to a 
width r/2 = 0.75 ± 0.2 GHz of the electronic quan¬ 
tum dot states. We assume that the main contributions 
to the broadening of the electronic states are the dou¬ 
ble dot decoherence and the coupling to the leads. In 
the emission measurements the decoherence lies between 
0.5 and 1 GHz while the coupling rate to the leads is 
Fl « Fij « 1 GHz. Because both of these contributions 
are of the same order of magnitude we cannot distinguish 
their individual effects. Investigating the effect of deco¬ 
herence on the emission linewidth would, however, be an 
interesting subject of further study. 
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